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Abstract. Even-even nuclei in the A ^ 100 mass region are investigated within 
the framework of the interacting boson model- 1 (IBM-1). The study includes energy 
spectra and electric quadrupole transition properties of zirconium, molybdenum, 
ruthenium and palladium isotopes with neutron number A^ > 52. A global 
y . parametrization of the IBM-1 Hamiltonian is found leading to a description of about 

pH I 300 collective levels in 30 nuclei with a root-mean-square deviation from the observed 

level energies of 120 keV. The importance of the (^5/2 subshell closure at neutron 
number A^ = 56 is pointed out. The geometric character of the nuclei can be visualized 
by plotting the potential energy surface V{f3, 7) obtained from the IBM-1 Hamiltonian 
in the classical limit. The parametrization established on the basis of known elements 
is used to predict properties of the unknown, neutron-rich isotopes ^°^Zr, ^^^Mo, ^^^Ru 
^ ■ and i22pd. 
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1. Introduction 

The production and use of radioactive beams is a rapidly developing new field in nuclear 
physics. Pioneering experiments are taking place, dedicated facilities are commissioned 
and new ones are planned. The primary aim of this research activity is the study of 
so-called exotic nuclei, that is, short-lived nuclei far removed from the line of stability 
because of an unusual ratio of numbers of neutrons A^ over protons Z. Parallel with the 
rapid developments of experimental techniques to observe exotic nuclei, there has been 
a surge in the interest in their theoretical description. 

The main motivation for these studies is the possibility that wide variations in 
N/Z may have a profound influence on shell structure, which constitutes the basis of 
our understanding of nuclei. There exist currently two standard theoretical approaches 
for the description of exotic nuclei, namely energy density functional theory [1] and the 
nuclear shell model |2]. In this paper we explore the use of a third approach, namely 
the interacting boson model (IBM) of Arima and lachello [Sj HI [5], to predict properties 
of exotic nuclei. 

The IBM is a semi-microscopic model of nuclei, positioned intermediately between 
single-particle and collective models. The model contains a vibrational and a rotational 
limit (as well as one which can be considered as intermediate) — which connects well 
with the phenomenology of nuclei — and can be brought into relation with the nuclear 
shell model. In the original version of the IBM, applicable to even-even nuclei, the 
basic building blocks are s and d bosons with angular momenta L = and L = 2, 
respectively. The s and d bosons can be interpreted as correlated Cooper pairs formed 
by two nucleons in the valence shell coupled to angular momenta L = and L = 2. This 
interpretation constitutes the basis of the connection between the IBM and the nuclear 
shell model, and leads in a natural way to a version of the model where a distinction is 
made between neutron and proton bosons, the so-called IBM-2 [6]. 

If one restricts oneself to configurations that are symmetric in the neutron and 
proton bosons, the simplest version of the model arises, known as IBM-1. The standard 
IBM-1 Hamiltonian has 1 single-particle energy and 5 two-body interactions if only 
relative (excitation) energies are considered, increasing to 2 + 7 if absolute binding 
energies are included in the analysis. These coefficients enter as free parameters in 
the model and must be adjusted to the data. However, because of the connection 
with the nuclear shell model, the dependence of the energies and interactions on the 
number of valence neutrons and protons is known qualitatively [7]. The long-term aim 
of this work is to develop a systematic IBM-1 parametrization (that is, to determine 
the functional dependence of the various coefficients on the number of valence neutrons 
and protons) for entire regions of nuclei, and to do so on the basis of a fit to the data, 
guided by microscopic arguments. Once the appropriate functional dependence has been 
determined, it will be possible to make spectroscopic predictions by extrapolation to 
unknown regions of nuclei. The predictive power of this procedure will depend crucially 
on the counting of the valence nucleons and hence on the shell structure of nuclei and 
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its evolution in exotic regions. 

The idea to use the IBM to extrapolate from known to unknown regions of the 
nuclear chart is not new. In fact, in the late 1970s and early 1980s, long before any 
talk of exotic nuclei, this technique was applied on many isotope series in the context of 
the IBM-2. For a review of the literature, see lachello and Arima [8]. Occasionally also 
isotone series have been considered |9] and recently also a combination of isobars and 
isotones [lOj. However, to our knowledge a comprehensive and systematic description 
with the IBM of entire regions of nuclei has not yet been tried. The closest related 
study is due to Casten et al. [11] who elaborated a multi-nucleus IBM-1 parametrization 
based on the N^Nt^ scheme [12]. Our intention is to develop this idea further with 
more refined Hamiltonians that better reproduce observed nuclear properties. We also 
mention that our investigation has similar aims as those of the recent work of the Tokyo 
group [13], {i.e., the use of the IBM in the description of exotic nuclei); in the latter 
studies the connection of the IBM with mean-field models is exploited to arrive at a 
suitable parametrization while here the long-term aim is to draw inspiration from its 
connection with the shell model. 

In this paper we present a first calculation of this type and apply the IBM-1 to an 
entire region of the nuclear mass chart, namely the even-even Zr, Mo, Ru and Pd nuclei 
with A ~ 100. Section [2] gives a description of the parametrization of the Hamiltonian 
used. The application to the A ~ 100 region is discussed in section [3l In the final 
section H] the conclusions and outlook are presented. Preliminary versions of this work 
were presented in references [HI [15] . 

2. Parametrization of the Hamiltonian 

The Hamiltonian used in this paper is of the form 

H = end + kQ ■ Q + k' L ■ L + k" P+ ■ P^ + \hl (1) 

where fid, -P+, L and Q are the boson- number, pairing, angular momentum and 
quadrupole operators, defined as 

hd = V5[S X rf]f , 

P+ = [st X st + V5d^ X d^S'\ P- = {P+y , 

L^ = VlO[d^ X J]«, 

Q^ = [d^xs + s^ X rf]f + xM^ X J]f . (2) 

Equation ([1]) defines an IBM-1 Hamiltonian in terms of the six parameters e, k, k', k", 
A and x- If one is interested in the properties of a single nucleus, this is a possible 
form of the most general Hamiltonian with the advantage that the first four terms have 
been used extensively in phenomenological fits to nuclear spectra [Sj. The last term 
hd gives rise to a so-called 'r-compression' which increases the moment of inertia with 
increasing angular momentum (or with increasing rf-boson seniority r) [16] . In addition 
to measured energy spectra, also electric quadrupole transition rates are considered in 
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the fit. These are calculated in IBM-1 using the E2 operator 7)j(E2) = e\^Q^ where Cb is 
the effective charge of a boson and with a quadrupole operator which is consistent with 
the one used in the Hamiltonian, following the consistent-Q formalism (CQF) [T7]. 

Although reasonable results are obtained with a constant Hamiltonian, they are 
considerably improved if parameters are allowed to vary with the number of nucleons 
in the valence shell, as suggested by the shell- model interpretation of the IBM-1. We 
propose here a dependence of the generic form 

X = Y.^^AU)\h)\ (3) 

where x is a parameter of the Hamiltonian, that is, x = e, k, k', k" or A, and fp is 
the fractional filling of the neutron (p = v) or the proton (p = vr) valence shell, that 
is, fp = Up/Qp with Up the number of valence neutrons or protons and flp the size 
of the corresponding valence shell. Note that in the parametrization ([3]) neutrons and 
protons are always counted as particles and never as holes. While microscopic arguments 
suggest a dependence on Fp = fip/VLp rather than on fp (where fip is the number 
of particles or holes, whichever is smaller counted from the nearest closed shell), the 
latter parametrization has the drawback of yielding spectra that are symmetric around 
mid-shell. Since this particle-hole symmetry is broken in the nuclei considered here, 
especially in the Pd isotopes, we prefer the parametrization (j3]). 

3. Application to the A ~ 100 region 

In the present application to the A ~ 100 region a linear dependence on the fractional 
fillings is assumed for the parameters k, k' and A while e and k" are kept constant. 
First, energy levels in the known nuclei ^^'^-^^^Zi, ^"^"^^^Mo, ^"^"^^^Ru and i02-i20pj ^^^ 
selected and fitted with the IBM-1 Hamiltonian ([1]) by minimizing the root-mean-square 
(rms) deviation. In this first step an initial choice is made for x, and e, k, k', k" and A 
are fitted to the energy spectra of all nuclei considered in the fit with no dependence on 
the number of neutrons or protons. Next, the linear dependence on the fractional filling 
of the parameters k, h' and A is determined again by minimizing the rms deviation. 
The minimization in this case involves 11 parameters, three for each coefficient in the 
Hamiltonian with a fractional-filling dependence together with e and /t". Once the wave 
functions of all states have been determined in this unified fit, the E2 transition rates 
are fitted by adjusting the value of the boson effective charge Cb but keeping to the 
same value of x- Since i?(E2) properties depend rather strongly on x, it is at this point 
that it can established whether its value is appropriate or not. The entire procedure 
is repeated for a new value of Xi until a reasonable compromise is found between the 
results of the energy and E2 fits. 

This procedure is illustrated in table [U The rms deviations in the energies and 
B{Ei2) values are quantified with following definitions: 
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Table 1. Root-mean-square deviations A in the energies and B(E2) values as a 
function of the quadrupole parameter x. 



X -0.05 -0.1 -0.2 -0.3 


-0.4 -0.6 -0.8 


-1.0 




A{E) (keV) 126 122 120 122 

Ai(E2) (e^b^) 0.26 0.13 0.10 0.11 

A2(E2) 3.44 2.29 2.02 2.12 


125 129 131 
0.14 0.17 0.15 
2.34 2.54 2.20 


133 
0.16 
2.21 




Table 2. Parameters of the Hamihonian ([T]) 


m units of keV. 






e kqo Kio kqi Kqq Kj^q Kqi k, Aqo 


Aio 


Aoi 


695 -35.8 -72.4 47.9 18.3 12.8 - 


15.8 1.4 -115.9 - 


-166.1 


161.2 



Ai(E2) 



\ 



iV, 



E2 



Y: (5(E2)L - B{E2)l^f 



A2(E2) 



^ ( B{E2)l^-B{E2)\A 
\N^,^\ 5(E2)^, j ' ^ ^ 

where the sums are over the available data points, A'^^; and A''e2 in number, respectively. 
Note that the first rms deviation for the E2 data, Ai(E2), is rather insensitive to small E2 
transitions while A2(E2) probably assigns too much weight to them. The rms deviations 
vary rather weakly except for small values of x, but both energies and E2 transitions are 
optimized for x ~ —0.20. This value of x ci-nd the parameters given in table [2] determine 
all spectra shown below. 

In the above procedure care has been taken to undertake a gradual release of the 
parameters, instead of immediately attempting a full 11-parameter fit. This method 
has proven to be numerically stable but it offers no guarantee that the global minimum 
of the 11-dimensional parameter space will be reached, instead of a local one. Finding 
the global minimum is all the more problematic since some of the parameters, especially 
those that account for the fractional-filling dependence, offer little intuitive insight and 
are thus entirely unknown a priori. The character of the minimum corresponding to 
the parameters of table [2] has been checked by choosing different starting values for the 
Xoi and Xio parameters, leading always to final parameter sets that are close to the one 
shown, within numerical errors. While this is not a proof of the global character of the 
minimum, it is a good indication of it. 

The experimental and calculated energy spectra are compared in figures [1] and 
[21 The experimental energies are taken from the National Nuclear Data Center [18] 
with additional information concerning the nuclei ^°*Pd and ii2,ii8,i20pj^ from the 
references p^ [20]. Some of the spin assignments are uncertain (brackets in NNDC); 
these are not shown on the figure to avoid overcrowding. In the even-even Mo, Ru and 
Pd isotopes with N > 52, all known levels of the ground-state, quasi-7 and quasi-/? 
bands with angular momentum up to J = 10 are considered in the fit while for Zr only 
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Figure 1. The experimental (symbols) and calculated (lines) energy spectra for the 
Zr isotopes with 52 < TV < 66 (left) and the Mo isotopes with 52 < iV < 70 (right): 
(a) ground-state, (b) quasi-7 and (c) quasi-/3 bands. Only levels in Zr isotopes with 
A^ > 60 are fitted. All experimental levels shown for Mo are included in the fit. 



IBM for Ar^lOO nuclei 




54 56 58 60 62 64 66 68 70 72 74 N 



56 58 60 62 64 66 68 70 72 74 76 N 



>3 



!h Z 

c 
W 




Ru(b) 



10+ 




9+ 


,_^ 




^3 


8+ 


^ 


7+ 






>^ 


6+ 


2^2 


5+ 


c 


w 


4+ 




3+ 


1 



2+ 



-+- 



-+- 



54 56 58 60 62 64 66 68 70 72 74 N 



>3 



>, 

£?2 

u 

c 
W 




6+ 

4+ 


£?2 
c 

w 


2+ 




0+ 


1 



OLh ^- 



>3 




54 56 58 60 62 64 66 68 70 72 74 N 



56 58 60 62 64 66 68 70 72 74 76 N 



Figure 2. The experimental (symbols) and calculated (lines) energy spectra for the 
Ru isotopes with 54 < TV < 72 (left) and for the Pd isotopes with 56 < iV < 76 (right): 
(a) ground-state, (b) quasi-7 and (c) quasi-/? bands. All experimental levels shown are 
included in the fit. 
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Table 3. Root-mean-square deviations A(_E) when Zr isotopes with neutron number 
N > No are inchided in the fit. 



A^o 


60 


58 


56 


54 


52 


A{E) (keV) 


120 


158 


214 


215 


217 



isotopes with A^ > 60 are included. Figure [T] nevertheless also shows the experimental 
levels for ^"^-^^Zt {not included in the fit) merely to indicate that these strongly deviate 
from what is predicted from systematics. To give an illustration of the deviations, the 2f 
level in ^^Zr is observed at 1751 keV while it is calculated at 444 keV. These differences 
are indicative of the closure of the 1^5/2 shell at A^ = 56, resulting in a nearly doubly- 
magic behaviour of ^^Zr. These considerations can be made quantitative by including 
the ground-state bands of some of the lighter Zr isotopes in the energy fit. As illustrated 
in table [3] even the inclusion of just one isotope, ^'^Zr, yields a dramatic increase in the 
rms deviation of the fit to the entire region. Effects of 1^5/2 shell persist in ^^''Zr which 
exhibits a low-lying K^ = O"*" band at 331 keV. This clearly cannot be the collective 
first-excited K^ = O"*" band of the IBM-1 and therefore the band is not included in the 
fit. The next-excited K'" = 0+ band at 829 keV is possibly the quasi-/? band of the 
IBM-1, which is calculated at 682 keV. 

We emphasize that the IBM-1 parametrization ([3]) depends on the counting of the 
number of valence nucleons and hence on the definition of shell closures and magic 
numbers. Since A^ = 56 is not included as a magic number, it is therefore no surprise 
that the structural features associated with this shell closure are absent from the present 
calculations. The inclusion of all Zr isotopes in a successful fit would necessarily require 
a refined parametrization that takes account of the A^ = 56 closure. However, even at 
the present level of approximation, we believe to have illustrated with this example that 
the unified fits of the type used here are sensitive to nuclear shell effects. This feature 
could be of use to detect precursor effects of shell closures in extrapolations toward 
unknown regions of the nuclear chart. 

In the present application excluding Zr isotopes with A^ < 60, we arrive at a 
description of 335 collective levels in 31 nuclei with a rms deviation from the observed 
level energies of 120 keV. The most notable systematic discrepancies are found for (i) 
the spectra of A^ ~ 56 Mo isotopes due to some remnant effect of the A^ = 56 closure, (ii) 
the energy of the 7 band in ^°^Zr which is predicted lower than it is observed, and (iii) 
some mid-shell nuclei like ^°°Ru and ^°^Pd which behave less collectively than predicted 
by IBM-1 systematics. 

On the basis of the parametrization thus established it is now possible to make 
predictions for more exotic nuclei for which up to now no spectroscopic information is 
available. This is illustrated in figures [1] and [2] where systematically the predictions for 
the next (unknown) even-even isotopes are shown, that is, for ^°^Zr, ^^^Mo, ^^^Ru and 
^^^Pd. Clearly, predictions further away from stability are possible as well but become 
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Figure 3. The experimental (full symbols) and calculated (open symbols) i?(E2; l\ 
0^) values in the Zr, Mo, Ru and Pd isotopes. 



increasingly uncertain. 

Once the wave functions of all states have been determined in this unified fit, one 
can also calculate the E2 transition rates with an overall effective charge Cb = 0.097 eb, 
adjusted to reproduce the iX — t- {s\ and ^ — )■ iX transitions. The results for the 
B{EI1\2X — !■ 0]^) and 5(E2;4^ — )• 2X) values are shown in figures |3] and HI respectively. 
These are a measure of the collective behaviour of the nucleus and our calculation 
generally follows the observed trends. A glaring exception to this overall agreement is 
the 2X — )■ O^*" transition in the Zr isotopes with A^ < 60 where the observed i?(E2) value 
is only fraction from what is expected from systematics. This again is due to the neglect 
in the calculation of the 1^5/2 closure at A^ = 56. There is also a mid-shell mismatch for 
the Pd isotopes which is calculated at A^ = 66 but observed at A^ = 64. Figure O shows 
the i?(E2;0^ — ?■ iX) values. In the vibrational U(5) limit this transition corresponds 
to one from a two-phonon to a one-phonon state which is allowed while in the more 
deformed limits S0(6) and SU(3) this i?(E2) value is much smaller. This argument 
gives a qualitative understanding of the theoretical results shown in figure O 

As a further illustration of the kind of accuracy obtained for the i?(E2) values in 
this overall fit, we show in table H] the results for three nuclei, ^"''Ru, ^°^Pd and ^"^^Pd, 
for which extensive E2 data is known. 

The geometric character of all nuclei can be visualized by plotting the potential 
energy surface y(/3,7) obtained from the IBM-1 Hamiltonian in the classical limit 
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Figure 4. The experimental (full symbols) and calculated (open symbols) -B(E2; 4^ 
2^) values in the Zr, Mo, Ru and Pd isotopes. 
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Figure 5. The experimental (full symbols) and calculated (open symbols) -B(E2; O2 
2^) values in the Zr, Mo, Ru and Pd isotopes. 
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Figure 6. Potential energy surfaces for all nuclei (first part). 
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Figure 7. Potential energy surfaces for all nuclei (second part). The grey background 
indicates a potential surface for a predicted nucleus. 
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Table 4. Experimental and calculated -B(E2; J^ -^ J[) values in units of 10^ e^fm"' 
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.[23]. The classical limit of the IBM-1 Hamiltonian ([T]) is 



"^^"^ ~ '^ ^^''' ' aS/3^ + ai?/3^ cos 87 + 4?/?^) , (5) 



:i + /32)' v-oo 
where the coefficients a^" are given by 

aff = Sk, ajo^ = e + {l + x^)i^ + QK' + A, afij = -k", 

ag^=4«:-^«:", 4? = -4^^x'^, afj = ^x^t^ +\'^" + ^- (6) 

Since the parameters in this potential depend on the fractional ffilings of the valence 
shells, we arrive at a geometry for each nucleus which changes with neutron and proton 
numbers and which has been obtained in an unbiased and straightforward way from the 
observed excitation spectra. This is illustrated in figures [6] and [71 where the potentials 
are shown for all nuclei included in the fit as well as for some which are predicted. For the 
sake of a correct interpretation of these potential surfaces, we recall that 7 is identical 
to the corresponding parameter in the Bohr-Mottelson geometric model [24] while /3 is 
only proportional to it, the proportionality factor depending on the ratio of the valence 
over the total nucleon numbers [25]. The potential energy surface of the nucleus ^"^Ru 
has recently been calculated in variety of ways. Moller et al. [26] showed that deviations 
from axial symmetry can be important and, in the specific case of ^°^Ru, obtained a 
potential energy surface with triaxial minimum. With the simple parametrization of 
the IBM-1 Hamiltonian as proposed here which is limited to two-body interactions, no 
triaxial minimum can be obtained [27], a result confirmed by the surfaces shown in 
figures [6] and [71 It is also seen, however, that the specific potential for ^°*Ru is rather 
soft in 7. Triaxial minima can be obtained when cubic terms are added to the IBM-1 
Hamiltonian [2H1 [29]. In comparison with the results obtained in reference [29] for the 
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nuclei ^°^'^^°'^^^Ru, we observe that the potential surfaces shown here show less variation 
with neutron number since they are derived from a global fit to many nuclei while in 
the former study the fits were carried out nucleus by nucleus. Also shown in figures [6] 
and [7| are the potential energy surfaces for the predicted neutron-rich nuclei. As can be 
seen from the figures, their geometry varies from strongly deformed in the neutron-rich 
Zr isotopes to spherical in ^^^Pd. 

4. Conclusions 

We have suggested in this paper the use of IBM-1 for global calculations of large regions 
of the nuclear chart and have illustrated the idea with an application to the A ~ 100 
region. The present work calls for several features to be studied in more detail. The first 
concerns the technical issue of finding the set of parameters that corresponds to the global 
minimum of the rms deviation. The minimization algorithm that has been used here is 
a straightforward linearization procedure of the eigenvalue problem associated with the 
Hamiltonian ([T]) which has no trouble in finding a local minimum. We have found no 
dependence of the converged parameters on various choices of their initial values but it 
would be of interest to confirm the global character of the minimum using a different 
algorithm which samples large parameter regions. Secondly, the dependence of the 
parameters on the fractional fillings fp of the valence neutron and proton shells requires 
further study. We have proposed here a dependence on Up/VLp but other choices should 
be investigated as well, such as a dependence on Np/flp (where Np is the boson number) 
or on Casten's factor P = Nt,NT,/{N,^ + N.„) [I2]. Furthermore, a possible dependence on 
subshells must also be studied. A proper treatment of these effects and their dependence 
on valence nucleon number must draw inspiration from the microscopic foundation of 
the IBM. For example, in this study we have taken a constant value of x, the parameter 
appearing in the quadrupole operator. This is possibly a reasonable approximation 
for the restricted set of nuclei considered here but it will certainly be inadequate if 
large regions of the nuclear chart are fitted simultaneously. This inadequacy might be 
satisfactorily resolved by a proper estimate of this parameter based on the shell-model 
interpretation of the IBM. An alternative strategy is to adopt a functional dependence 
of X on valence-particle number suggested by microscopy, containing a few parameters 
that are fitted to the data. Finally, we point out that the type of calculation described 
here opens up the possibility for a simultaneous and global treatment of nuclear ground- 
state properties (such as masses and radii) in addition to the excited-state properties 
discussed in this paper. A possible strategy for merging the calculations of ground- 
and excited-state energies in the framework of the interacting boson model was outlined 
in reference ^U\ with preliminary results for even-even nuclei in the major shell with 
82 < A^ < 126 and 50 < Z < 82 but further investigations of this approach are required. 
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